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We present the theory of the effect of thermal fluctuations on commensurate p x p density wave 
ordering on the square lattice {p > 3, integer). For the case in which this order is lost by a second 
order transition, we argue that the adjacent state is generically an incommensurate striped state, 
with commensurate p-periodic long range order along one direction, and incommensurate quasi-long- 
range order along the orthogonal direction. We also present the routes by which the fully disordered 
high temperature state can be reached. For p = 4, and at special commensurate densities, the 4x4 
commensurate state can melt directly into the disordered state via a self-dual critical point with 
non-universal exponents. 



I. INTRODUCTION 

A variety of remarkable recent low temperature scan- 
ning tunneling microscopy (STM) observations have re- 
vealed periodic modulations in the local density of states 
of the cuprate superconductorsiiSiii. Notably, the ob- 
servations of Hanaguri et a& Ca2-a;Naa;Cu02Cl2 clearly 
show that the modulations have a commensurate period 
of 4 lattice spacings along both directions of the under- 
lying square lattice (4x4 ordering). Closely related, but 
not identical, modulations have been observed in higher 
temperature STlVfi, and in neutron scatterin g^'^' ^ ob- 
servations. The differences between the various exper- 
iments relate (i) to the period of the ordering, which 
is incommensurate in Ref. Hi and (ii) to whether the 
ordering extends along one or both of the x and y 
axes — neutron scattering experiments are more easily ex- 
plained by anisotropic ordering along one of the axes 
directionsiSiiiiiS. An important open question is whether 
these differences reflect a fundamental distinction in the 
underlying electronic structure of the cuprates, or they 
can be explained by the differing experimental parame- 
ters of temperature and carrier concentration. 

This paper will begin with the assumption that at very 
low T, for a small range of carrier densities, and in a suffi- 
ciently clean sample, the system has perfect 4x4 density 
wave order (and, more generally, p x p order with p > 3, 
integer). We will then present a general phenomenolog- 
ical theory of the effect of thermal fluctuations on such 
an ordered state, describing the phases that appear at 
higher T and for a wide range of carrier concentration. 
Because we focus exclusively on thermal fluctuations, our 
results are quite general and independent of the precise 
microscopic nature of the density wave ordering; in par- 
ticular, the ordering could be a site charge density wave, 
or a modulation in exchange or pairing energies due to 
valence bond solid order. Our results depend only on the 
fact that spin-singlet observables acquire a periodic mod- 
ulation at low T, and follow completely from the symme- 
try properties of such states. The thermal fluctuations 
of such an ordered state can, and will, be described by a 
purely classical theory of the density wave order param- 
eters. 



There is a actually large early literature on the melt- 
ing of a variety to two-dimensional solids on different 
substrates, and on the commensurate-incommensurate 
transition"d'^d^d^d'^d^^Qi^°. However, most of this 
work considered the case of the triangular solid on sub- 
strates with six-fold symmetry, or the commensurate- 
incommensurate transition on anisotropic solids. The 
particular case of interest here, a square lattice substrate 
and isotropic p x p ordering, appears not to have been 
considered previously. We will therefore present an ex- 
tension of this earlier theory to the situation appropriate 
for the cuprate superconductors. As we will show below, 
our results have a number of novel features not found in 
the cases studied earlier. 

We will begin in Section^ by deflning the order pa- 
rameters of the commensurate p x p solid on the square 
lattice. Symmetry considerations then allow us to ob- 
tain a phenomenological free energy density controlling 
thermal fluctuations of the order parameter, and a corre- 
sponding mean fleld phase diagram shown in Fig. ^ At 
this stage, this phase diagram contains 3 phases whose 
characteristics we summarize below: 

(A) Commensurate solid: This is the ground state 
with long-range density wave order with period p along 
both the X and y axes. 

(B) Incommensurate, floating solid: In mean fleld 
theory, this phase has density waves with the same in- 
commensurate period along both the x and y directions. 
This order is only quasi-long-range i.e. density wave cor- 
relations decay with a power-law at long distances. 

(C) Liquid: This is the high temperature phase in which 
all density wave correlations decay exponentially with 
distance. 

The remainder of the paper will consider fluctuation 
corrections to this mean field phase diagram. 

In mean field theory, the transition between phases 
A and B is first order, with a jump in the period 
of the density wave from a commensurate to an in- 
commensurate value. We know from previous work 
on anisotropic phasesi^^i^iiSiiS that such commensurate- 
incommensurate transitions are driven by the prolifera- 
tion of domain walls, and we consider the domain wall 
theory for such a transition in Section IIIII In our case 
there are two sets of domain walls, running predomi- 
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FIG. 1: Mean field phase diagram of the fi'ee energy T,s, in 
Eq. 12.41 . Ahhough all parameters in T,f, axe functions of 
both temperature and carrier concentration, the vertical axis 
above, a, is primarily a function of temperature. The hori- 
zontal axis is tuned by rj, which is mainly a measure of carrier 
concentration away from the commensurate value appropriate 
to the p X p ordered state. There are three stable mean field 
phases separated by a first order transition (thick line) , a sec- 
ond order transition (solid line), and a second-order critical 
point at s = 77 = 0. The phases A, B, C, are described in the 
text. For phases A and B we show a schematic of the density 
modulation 6p{r) for the special case p = 4, with the points 
representing the underlying lattice. Upon including fiuctu- 
ations, the transition between phases A and B can become 
second-order; in this situation the phase adjacent to phase A 
is a new phase Bs, the incommensurate striped phase. The 
phase Bs preempts a portion (or all) of phase B. See Fig|l| 
for the fiuctuation-corrected phase diagram of phases A and 
B. 



nantly in the x and y directions, and a key parame- 
ter will be the intersection energy, //, of these domain 
walls. If // < 0, then the mean field prediction of a first 
order transition to an isotropic floating solid is main- 
tained. However, for the case of a positive intersection 
energy (// > 0), we will demonstrate that there is a 
second order transition to an anisotropic, incommensu- 
rate "striped" state, which we denote B5. The B5 state 
has commensurate long range order with a period of p 
lattice along one axis, and incommensurate quasi-long- 
range order along the orthogonal axis. Because Bg has 
long-range order only along one direction, it may also be 
labelled a 'smectic'iSLS^; however the presence of the un- 
derlying lattice makes its fluctuations quite different from 
a conventional smectic liquid crystal. The striped solid 
B5 will replace a portion of phase B, so that Bg covers 
at least the entire second-order phase boundary to phase 
A (see Fig 0J . The remainder of phase B could remain 
an isotropic floating solid, or it could be anisotropic in 
its entirety. 

Section IIVI will consider the transitions from the in- 
commensurate states in phase B to the disordered liquid 
C. Two distinct scenarios are possible, depending upon 



the whether the transition takes place from an isotropic 
floating solid B, or from the anisotropic solid (B5) which 
has incommensurate order along only one axis. In ci- 
ther case, the transition is driven by the unbinding of 
dislocationg-i?:!-'^^, and the complete theory for such tran- 
sitions will be presented. 

Finally, Section IVI will present the theory of the direct 
transition from the commensurate p x p solid A to the 
liquid state C. From Fig. it appears that such a tran- 
sition is only possible at special commensurate values of 
the density. For the physically relevant case of p = 4 we 
will find that a direct second order transition is possible. 
It is described by a self-dual theory with continuously 
varying exponents, and is a generalization of the theory 
found in Ref. for the XY model. 



II. ORDER PARAMETERS AND MEAN FIELD 
THEORY 

As noted in Section all phases and transitions ex- 
amined here are associated with the order of a generic 
'density', which could be any observable invariant under 
spin rotations and time reversal. We represent this den- 
sity by 6p{r). At sufficiently low temperatures, in weak 
disorder, and for some range of carrier concentration, we 
assume that this density prefers to order with a period 
of p lattice spacings (p = 4 is the case of interest for the 
cuprates). We can therefore write 

Sp{r) = Re [$^6'^==-'"] + Re [^yc'^y"] , (2.1) 

where = (27r/o)(l/p, 0), Ky = (27r/a)(0, 1/p), and 
^x.y are complex order parameters which vary slowly the 
on the scale of a lattice spacing. 

We now want to write down the most general free en- 
ergy for consistent with the symmetries of the un- 
derlying lattice. Among these are T^^y which translate 
by a lattice spacing in the x, y directions, and Ix,y which 
reflect the x, y axes. These operations lead to 
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We will also assume the symmetry of rotations by 90 
degrees, R, under which 

R:^.^% ; %^K- (2.3) 

This symmetry is absent in some of the cuprate com- 
pounds (most notably, in YBCO), and it is not difficult 
to extend our considerations to include this case. 

We can now write down the most general free energy 
density, expanded in powers of ^x.y and its gradients, 
consistent with the symmetries in Eqs. (|2.2|l and H2.3|) 
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(see also Refs. |2J and|25D. This is 



+ C2 (la,,*,!' + 

i;|$^n$y|2 + u;$P+c.c. + w$P+c.c. . (2.4) 



Here s is a parameter which we wih tune to drive the 
transition; it is assumed to contain the primary depen- 
dence on temperature, and wiU also have some depen- 
dence on carrier density. 

The term proportional to rj is allowed by the symme- 
tries. It indicates that at sufficiently high temperature, 
the density correlations are generically incommensurate. 
This incommensurability is a consequence of the values 
of the domain wall energies between the p commensu- 
rate ordered statesi8i2£: a domain wall between state 1 
and state 2 (say) will generally have a different energy 
than a domain wall between state 1 and state p. In other 
words, if we represent the p states along one direction by 
a p state clock model, then the interactions of the clock 
model are 'chiral'ii^ 

The sign of the parameter v implies a preference for 
either isotropic order {v < prefers both {^x) and ($y) 
non-zero), or anisotropic order (v > prefers only one of 
{^x) or ^y non-zero). We will assume throughout this 
paper that < 0, so that {^x) and ($j,) are both non-zero 
at sufficiently low temperatures. Nevertheless, we will 
find that thermal fluctuations can induce an anisotropic 
striped state over an intermediate temperature range. 

The complex parameter w accounts for the commensu- 
rability lock-in energy. It ensures that at sufhciently low 
temperatures (specifically, for s sufficiently negative), the 
order has a commensurate period of p lattice spacings. 
For V < (assumed throughout), the ordering will be 
px p. The phase of w controls the phase of this ordering 
('site-centered' or 'bond-centered'). None of our results 
will be sensitive to this phase, and will apply equally to 
all of them. 



A. Mean field theory 

Now we present the results of a simple mean field min- 
imization of leading to the phase diagram in Fig. ^ 
As described in the introduction we find three distinct 
phases. The low temperature ground state, labelled by 
A has long rage density wave order with period p char- 
acterized by 



where 9{p) = (2n -I- l)7r/p if the phase is bond-centered 
and 9{p) ~ 2mT /p if the phase is site-centered with n — 
0,l,...,p— 1. The inset in Fig. QA shows the locked in 
solid for the special case of a bond-centered phase with 
p = A where the positions of the underlying atoms are 
indicated by the light gray circles. This commensurate 
phase has a first order transition described by the line 
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\/2u + V — \/2u + V — pw J 4Ci 



(2.6) 



to B which has incommensurate or floating density wave 
order characterized by 
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The magnitude of both $^ and $y are equal and they 
have the same incommensurate period along perpendic- 
ular directions. For the special case of p = 4 the floating 
state is shown as an inset in Fig.QlB. The floating solid 
can melt via a second order phase transition defined by 
the line 
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to a fully disordered state C with ^ = 0. Finally, 
there is a tri-critical point for s = rj = where the com- 
mensurate solid A can melt directly to the liquid phase 
C. 



III. DOMAIN WALL THEORY OF THE 
COMMENSURATE-INCOMMENSURATE 
TRANSITION 

In Section ing we found a mean-field first order tran- 
sition between the commensurate p x p solid A and the 
incommensurate, isotropic floating solid B. This section 
will examine fluctuations near this transition more care- 
fully. Wc will find that the transition can actually be 
second order under suitable conditions, and the second- 
order order transition is to an anisotropic, striped state 
with commensurate p period order along one direction, 
and incommensurate order along the orthogonal direc- 
tion. 

For our study of the initial melting of the commensu- 
rate ordered state, we will assume that dislocations can 
be ignored. The effect of dislocations will be considered 
in Sect ion Hvl and we will then verify the self-consistency 
of this assumption. Instead, the primary actors will be 
domain walls between the commensurate states, as is 
also the case in previous thcories^^ of the commensurate- 
incommensurate transition in anisotropic systems. 

In the absence of dislocations, we can focus on globally 



defined single-valued angular variables 0x.y with 



^y cx e*' 



(3.1) 
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(2.5) 



Note that the values of 6x.y span over all real numbers, 
and not just modulo 2tt. However, periodic boundary 
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conditions need only be satisfied modulo 1-k. In the naive 
continuum limit, the free energy for d^ y can be expanded 
in powers of the local "strains" Vr^^.y: 



solution is a soliton with equation 



+ K4(d^9,){dy9y)-Tj[d^e^+dyey] 

— h [cos{p9x) + coa{p9y)] — . 



(3.2) 



All terms above are, in principle, obtained from those in 
Eq. H2.4|l . but now we have retained terms up to second 
order in spatial gradients. Now, the incommensuration 
is induced by the total derivative terms proportional to 
rj: these are non-zero because the angular fields can ac- 
cumulate an integer multiple of 27r even under periodic 
boundary conditions. Similarly, we do not have the free- 
dom to integrate by parts, and so combine the terms pro- 
portional to Kg, and K4. The commensurability energy 
is now imposed by the p-fold field hp. 

We will begin in Section IlII Al bv describing the mean- 
field structure of the domain wall (or 'soliton') excita- 
tions of J-g- Then, in Section [111 ij| we will present the 
theory of the commensurate-incommensurate transition 
driven by the proliferation of these domain walls. Sec- 
tion nil CI will address the nature of density fluctuations 
within the incommensurate phases: we will estimate the 
elastic constants of these phases using the domain wall 
theory of Section IIII Bl 



A. Energetics of domain walls 

Consider starting at low temperatures from a fully or- 
dered two-dimensional crystal. In this state 6^ — 27rn/p, 
9y = 27rn'/p, where n, n' are integers, at all points in 
space. We are now interested in the deviations from this 
perfectly ordered state as measured by the continuum 
free energy J^g in Eq. (|3.2I) . 

The simplest deviation is a single line domain wall in 
which 0^ increases by 27r/p across the domain wall. For 
r] > (which we assume, without loss of generality), such 
a domain wall will preferentially run in the y direction 
i.e. the domain wall has 9y constant, and 9x a function 
of X only. The x dependence of Ox can be determined by 
the sine- Gordon saddle point equation 



Kid%{x) ^ phsm[pex{x)] 
subject to the boundary conditions 

0,(0) = ; e,{Lx) = — 

P 



(3.3) 



(3.4) 



c{x) = — tan ^ 
P 



^p^h/Ki{x-L^/2) 



(3.5) 



which allows us to determine the value of the domain wall 
free energy per unit length (excluding the contribution of 
the r] term in J-^i), which we represent by e, 



P 



(3.6) 



Similarly, there is a corresponding domain wall in 9y, 
with identical physical properties. 

Now we consider the interesting case with domain walls 
running in both the x and y directions. These walls will 
intersect, and we are interested in the nature of the in- 
tersection, and of the intersection free energy //. We can 
focus in on a single domain wall crossing by imposing the 
boundary conditions of Eq. (|3.4(l in both the x and y di- 
rections. For the special case of perfectly straight domain 
walls, 9x{x,y) = 9x{x) and 9y{x,y) = 9y{y), the only in- 
teraction term in Eq. H3.2(l can be integrated exactly to 
give 



(3.7) 



where is the size of the system in the x-direction. The 



indicating that the sign of the domain wall intersection 
energy is equal to the sign of K4. When the domain walls 
are not straight, the interaction energy can be found by 
numerically minimizing J-g and evaluating the interaction 
terms at the ground state field configurations as seen in 
Fig. 121 For the following discussion we refer to Fig.|21and 
it is assumed that all stiffnesses are measured in units of 
Ki. For positive values of K4 and K2 = 1.0 (panels a 
through c) we observe straight domain walls for K3 < K2 
and only observe deviation when ~ K2. In panel d. 
Kg ~ K2 < K4 and we find that the walls wander quite 
significantly. Panels a through d all have K4 > and in 
agreement with our earlier discussion, // > in all four 
configurations. Panel e has no interaction terms {Kg — 
K4 — 0) and consequently // — 0. For negative values 
of K4 we find that the domain wall intersection energy 
changes to a negative value. For IKg] < K2 the walls 
remain straight (panels f and g) but as Kg ~ K4 ~ —K2 
the intersection energy becomes large and negative and 
the system attempts to extend the domain wall overlap 
over a finite region (panel h) . 

Although we have only presented eight distinct config- 
urations here, all of the Kj parameter space was exam- 
ined. When the magnitude of Kg and K4 are small with 
respect to K2 and Ki we always find configurations with 
domain walls crossing at right angles where the domain 
wall intersection energy can be calculated using Eq. (|3.7|l . 



B. Proliferation of domain walls 

Now we imagine increasing the parameter rj so that 
the total free energy per unit length of a domain wall is 



5 



a. 


I 


— 

K^ = 1.0 
K3 = 0.2 
K, = 0.8 


h 

fj= 1.97 



g- 



c. 

i 




K^ = 1.0 
K3 = 0.2 
K, = 0.2 


f, = 0.49 



e. 




= 1.0 
K3 = 0.0 
K, = 0.0 


f, = 0.0 



-- 1.0 
= 0.2 
= -0.8 



f, = -1.97 



D. 




K, = 1.0 


r 


= 0.8 




K, = 0.8 

4 


fj= 1.59 




f. 




K, = 1.0 
= -0.2 
K , = -0.2 


f, = -0.49 




/ h P 

Ik,' 



FIG. 2: Eight possible configurations of the free energy den- 
sity for domain wall crossings depending on the numerical 
values of K2, K3 and K4, with p — 4. The grey scale plots 
the size of the local free energy density and all stiffnesses and 
energies are in units of Ki. Panels a. through h. show de- 
creasing domain wall intersection energy //, and for straight 
walls, the numerically computed value can be compared with 
(^V4)^4. 




FIG. 3: A net of wandering domain walls constructed using 
random walkers on a lattice with both hard-core repulsion 
and restricted phase space. If the domain walls are separated 
on average by a distance £x in the x direction and iy in the 
y direction, then we observe collisions between domain walls 
running in the same direction with separation on the order of 
el or el (e^ = ey here). 



(Fig. 13). This can be computed using the elegant free 
fermion mapping of Pokrovsky and Talapo"\ii^, which ap- 
plies here (essentially unchanged) separately to the do- 
main walls in each direction. Briefly, the argument runs 
as follows. Let Ux{y) represent the x co-ordinate of a do- 
main wall in Ox- Any y dependence in increases the 
total length of the domain wall, and so leads to a free 
energy cost 



1 



dy 



dux 
dy 



(3.9) 



eventually negative. In such a situation we expect a pro- 
liferation of domain walls, leading to the appearance of a 
floating solid with incommensurate density correlations. 
This section will discuss the theory of such a transition. 

Let the incommensurate state have domain walls in 9x 
with an average spacing tx , and domain walls in 9y with 
an average spacing 9y. From the energy per unit length 
of these domain walls, and their intersection energy, com- 
puted in Section IIII Al there is a clearly a contribution 
to the free energy per unit area, Td, given by (see Fig.j^J 
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(3.8) 



However, in addition to the simple energetic contri- 
butions in Eq. (|3.8|l . we also have to consider the en- 
tropic contribution of the wandering of the domain walls 



Now the partition function at a temperature T of such 
domain walls can be mapped onto that for free fermions 
with a density l/£x and mass e/T. From the ground state 
energy of such a free fermion system, we then obtain the 
contribution of the wandering of the domain walls to the 
free energy density 



6e 



1 



(3.10) 



We are now faced with the simple problem of minimiz- 
ing the free energy density 



^d(^x,ty) — J^d \^x,(-y) + ^d \^x,^y) 



(3.11) 



and £y to determine the nature of the 



as a function of 
commensurate-incommensurate transition. This simple 
calculation turns out to have some interesting structure 
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FIG. 4: Phase diagram of the commensurate-incommensurate 
transition obtained by the minimization of Eq. 13.1111 over the 
values of the mean domain wall spacing and . 



which we will now describe. The results of such a min- 
imization are summarized in Fig ^ The nature of the 
phases appearing here depends upon the sign of the do- 
main wall intersection energy //. 

For negative // corresponding to < in J-'g, as 
we increase the value of 77 the mean field prediction is 
preserved, and the ground state with long range charge 
density wave order melts via a first order transition to an 
isotropic incommensurate floating solid which we labelled 
as B. This transition occurs when 
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That is, for 77 < rjj^B the system is commensurate with 
= = but for 7] > rjj^g there is a jump to the 
isotropic floating state B with domain wall densities 
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If the domain wall intersection energy is positive, im- 
plying that collisions are disfavored (// > 0, K4 > 0), 
then as we increase the value of 77 from zero the state 
remains commensurate with long range order and no do- 
main walls {l/£x = = 0) until 



pe 



(3.14) 



where there is a second order phase transition to a state 
with (say) 
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(3.15) 



which has domain walls running in either the x or y di- 
rection but not both. This anisotropic floating solid or 



incommensurate striped state which we label B5 has long 
range charge density wave order with period p in one di- 
rection and quasi long range order with power law decay 
in the perpendicular direction. Near the second-order 
transition between A and B5, insertion of Eq. H3.15|l into 

Eq. l|TTT|l shows that J^d {e - 2m-j/pf/'^. 

As 77 is increased further for // > 0, the phase Bg 
persists until 
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(3.16) 



where there is a first order transition to the isotropic 
incommensurate floating solid B with and £y still given 
by Eq. H3.13|l . In fact, the solution in Eq. 13.13|1 is a local 
minimum of the free energy for all 77 > pe/(27r). However, 
its free energy behaves like J^d ^ ~(e — '^t^vIpY' upon 
approaching phase A. Close enough to phase A, this free 
energy is always larger than the free energy for phase B5. 

Therefore, if the intersection energy for domain walls 
is positive, the melting of the commensurate solid A al- 
ways occurs via a second-order transition to the striped 
anisotropic state Bg. 



C. Renormalization of elastic constants in 
incommensurate phases 

Once the domain walls have proliferated in phases B 
and B5, the density correlations become incommensurate 
and appear at wavevectors which are shifted from the 
commensurate values in Eq. (|2.1|l . We are interested here 
in the long wavelength "spinwave" density fluctuations 
about this incommensurate state. 

First, let us consider the phase B. Rather than defining 
the spin-wave variables about the commensurate state as 
in Eq. (|3.1|) . we now need to look at fluctuations above 
an incommensurate ordered state. So now we write 



<^y cx e*2^2^/(P^«)e*®« 



(3.17) 



because the incommensurate ordering wavevectors are 
(27r/a)(l/p -ha/4,0) and (27r/a)(0, 1/p a/^y). We are 
interested in the effective action for these 9x,y on a coarse- 
grained scale much larger than £x,y This action can be 
deduced by the same symmetry arguments made in Sec- 
tion ^1 It is not difficult to see that this action in the 
floating phase B has the same structure as J-g in Eq. (|3.2|l . 
except that now 7/ = /i = 0: 
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^ [{dx0y)^ + {dye,)^] + K3{dx0y)(dy0x) 



K4idx9x)idy9y) 



(3.18) 
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For the phase B5, the 9y variables (say) are locked at 
their commensurate value, and so only the 9^ variables 
will contribute to the spin-wave theory: 



= I (fr 



(3.19) 



An important point is that the stiffnesses Ki^/^ in 
Eqs. H3.18|l and H3.19|l be strongly renormalized from the 
bare values in Eq. (|3.2(l . Following the analysis of Ref . ITtI 
we will now estimate their renormalized values K^_^ in 
terms of the parameters appearing in the domain wall 
free energies in Section IIII Bl (the superscript R will be 
used for clarity only in this subsection). 

First, imagine that we impose a uniform compressional 
strain dxOx on the floating solid. On the average, this will 
cause the domain walls to move closer to each other, and 
change the value of £x to £x — Six ■ Because 6x changes by 
2n/p across each domain wall, we conclude that dxOx = 
{2Tr/p)S£x/£x- The change in £x will cause a change in 
free energy that can be computed from Eq. I)3.11|l . and 
so we conclude 

P^-" ■^'^ (3.20) 



47r2 d£l 



where the derivative has to be computed at the equilib- 
rium values of ix, £y which minimize Eq. H3.11|l . This 
determines 



(3.21) 



in both phases B and B5. Similarly, we can apply a 
combined compressional strain in 9x and 9y and conclude 



K. 



47r2 dlxd£y 



(3.22) 



This is non-zero only in the floating phase B, and we 
obtain 



47r2 



phase B only. 



(3.23) 



To determine K2 , apply a small uniform shear strain 
to 9x by inducing a non-zero dy9x. This will move the 
domain walls in 9x such that the position of each domain 
wall obeys dyUx = {p£x/(2'K))dy9x. Inserting this into 
Eq. H3.9|l . we obtain 



p£x 

2tt 



47r2 ' 



(3.24) 



in both phases B and B5. 

Finally, to determine we need to apply shear 
strains to both 9x and 9y. Neither of them causes a net 
change in the density of domain walls, or in the number 
of intersections between the domain walls. Consequently 
there is no change to the free energy beyond that already 
accounted for by -ftT^*, and hence 



Ki^ = 0, 



(3.25) 



in both phases B and B5. 



IV. DISLOCATION MEDIATED MELTING OF 
FLOATING SOLIDS 

We will now consider the transition from the isotropic 
floating solid B and the striped floating solid B5 to the 
disordered liquid phase C. These transitions are driven 
by the unbinding of dislocations. 

With the parameters of the "spin- wave" theory at hand 
m Section HITCl we can address the energetics of the the 
dislocations. Wc will derive the effective action for the 
dislocations in Section FlV Al and then obtain the renor- 
malization group (RG) flow equations in Section [IV Bl 



A. Dislocation interactions 

For the most part, this section will be restricted to a 
discussion of dislocations in phase B. The simpler case 
of the anisotropic phase B5 is easily obtained by only 
including those terms arising from the fluctuations of 9x ■ 

Dislocations are simply 'vortices' in the angular fields 
9x y under which 



fdr.V.9x = 2nmxiry) 
j)dr-Vr9y — 2nmy{ry) 



(4.1) 



where mx,y are integers at the vortex (i.e. dislocation) 
site ry , and the integral is over a contour that encloses r„ . 
Each dislocation is therefore characterized by a doublet 
of integers (to2;,TOj,). 

To compute the interactions between these vortices, it 
is useful to define continuum vortex densities by 



TOi(r) = } m(r„)(5(r - r^,) 



(4.2) 



where i — x,y. Then, after transforming to momentum 
(k) space, the relation Eq. 14. 1|) can be written simply as 



2n 

ki9j = ki"dj + jj^iiki 



(4.3) 



where e is the antisymmetric tensor, and dj is an arbi- 
trary smooth angular field which has no vortices. We now 
insert Eq. (|4.3fl into the free energy in Eq. (|3.18() , and 
minimize with respect to . The result for the total free 
energy of the vortices is then 



47r2 D{kx,ky 



+ PyK.iKf - K!)) + |m,(k)|2 {klK,{Kl - Kj) 
+ klK2{Kl - KD) + mx{'k)my{~]^)kxky {K^Ks 



KlKi 



(4.4) 
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where 



D{K,ky) = K^K2{kl + klf 

+ klkl {{Ki - K^f - {K^ + K^f) . (4.5) 



Note that the values of if /L4 from Scction llll Cl are to be 
inserted into the expressions above; here, and henceforth, 
the superscript R has been dropped. The interaction 
between the vortices in now determined by transforming 
Eq. 14.411 back to real space. This takes the form 



)Vi\ 

v<iv' 

+ my{ry)my{vy>)V {\vy - 



(4.6) 



where 0(r) = arctan(?//2;) is the angle the vector r makes 
with the X axis. We will not need the explicit form of the 
interaction W in our subsequent analysis, and so we will 
not specify it explicitly; the interaction V is given by 



Appendix IbI note that the if 1-4 couplings appear upon 
taking the inverse of the matrix of couplings between the 
strains in Tb] consequently, the inverse expressions for 
the if 1-4 in terms of the if 1-4 have exactly the same 
structure as in Eq. H4.9|l . The parameters in the last line 
of Eq. H4.7II are given by 



in = 



gA(if.,^) 



(4.10) 



K1+K2 



V{cf>) 



K1K2 
47r2 



(ifl+if2)2-(if3 + if4)^ 



d(p\n{\ cos((/> - ((5)1) A(ifi, ip) 



where the calculation of Lg is described in Appendix ^ 
For ifs — K4 — 0, V{(p) can be evaluated in closed form: 



47r(ifiif2)i/2 



In 



if 



1/2 



if. 



1/2 



if 1 sin^ 



if 2 COS^ i 



(4.11) 



V{r,(t>) 



dip 



^ikr cos(0 — Ip) 



1 



Ki{Kl - Kl) cos2 ip + if 2 (if 1 - if|) sin^ p 
^ K1K2 + ((ifi - if2)2 - (ifg + KiY) sinVcos2 p 

2y2 .2. roo^ka^ 

1^ / ^ / — k — [cos(fc?' cos((?!) - p)) ~ 1] 



X A{K„ip) 



p2y2. 



■LoHr/a) + V{^) 



(4.7) 



B. Renormalization group flows 

With the knowledge of the interactions between the 
dislocations, the renormalization group equations can be 
derived by the methods already described in Refs. fisl 
and 113. We introduce a vortex fugacity, y = e~^''^'^ and 
examine the effect of integrating out pairs of dislocations 
in an expansion in powers of y. A simple and standard 
analysis shows that the flow equation for the fugacity is 



where we have inserted a soft cutoff using the lattice spac- 
ing a, and 



(4.12) 



A(if„(^) 



(4.8) 



if 2 cos^ p + if 1 sin p3 



KiK2+{ (if 1 - if2)2 - (ifa + if4)2 sin^ p cos^ p 



is a function of the new couplings ifi_4 defined by 



ifi = 



p2y2 



if2 



4^2 K^-Ki 



2rp2 



47r2 K(-K. 



if 2 = 



if 4 = 



p^T^ ifi 
4^2 Kf-Ki 



p2y2 



4^2 Kl-Kl 



(4.9) 



At this point these definitions of the if 1-4 may be viewed 
as arbitrary variables, but we will see later in Section Ivl 
and Appendix ^ that these are the couplings that ap- 
pear in a self-dual mapping of the theory J-'g. Also, from 



The renormalization of the ifi from the vortices can 
be coniputed by the method described in Appendix B 
of Ref. [l3- We compute the renormalization of the elas- 
tic constants, ifi by determining the contribution of the 
vortices to a two-point correlation of the strains. Such 
a procedure leads naturally to flow equations for the 'in- 
verse' or 'dual' if^ couplings, and we obtain 



difi 
~df 

dK2 

dl 

de 

dIU 
~df 



)cos2(Ae^(«/^ 



p y T I d(p sm" 0e 
'0 



2 rh.VW/T 



(4.13) 
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We can convert these into equations for the Ki and obtain 



DIRECT MELTING OF THE 
COMMENSURATE SOLID 



dKi 




d£ 


T 


dK2 




di 


T 






de 


T 


dKi 




d£ 


T 



(Kf + Kl) 



{Kl + Kl) 



)sin2 cAe^(*)/^ 



)cos2 Se^^^y^ 



K2K3 / d<?!)cos^0e^('^)/^ 
Jq 

K1K4 I d(?!)sin2(/)e^('^)/'^ 



(4.14) 



While complex in appearance, the flow equations 
Eq. (|4.12l) and (|4.14l) predict a melting transition of the 
floating phase B into the liquid phase C which is in a 
universality class closely related to that of the Kosterlitz- 
Thouless (KT) transition. The phase B is stable provided 
Lq > 87r/(p^T). Evaluating La at the values of the K^ 
determined in Section Fill CI we find that phase B is sta- 
ble everywhere for p — A towards an infinitesimal vortex 
fugacity for the elastic constants in the domain wall the- 
ory; this justifies the neglect of dislocations in our study 
of domain wall proliferation in Section IlIII The present 
equations Eq. H4.14|l show how phase B will ultimately 
become unstable to a vortex unbinding transition once 
the vortex fugacity becomes larger. The flow equations 
in Eq. (|4.14() imply singularities in the elastic constants 
just before the melting transition which can be computed 
as in Ref.liE Some of the universal amplitude ratios here 
will be different from the KT transition, but all exponents 
and critical singularities will be as in the KT transition. 

The flow equations for the melting of the anisotropic 
phase Bs can be easily obtained from Eq. H4.14|l simply 
by setting K3 = K4 = 0. In this case, the values of the 
elastic constants in Section IIII CI and Eq. H4.10|l imply 
that 



1/2 



TT 

T' 



(4.15) 



The equation for the vortex fugacity Eq. H4.12|l implies 
then that the vortex fugacity flows to zero as long asii 
p^ > 8, which is certainly the case for p > 3. Hence the 
floating phase B5 is also initially stable to dislocation 
unbinding. At higher temperatures, there can be a dis- 
location unbinding transition in the Kosterlitz-Thouless 
universality class. Note that in this sequence of transi- 
tions, even after such a transition has occurred from the 
phase B5, order has only been lost in x direction, and 
commensurate long-range order remains in 9y. So the 
system is stfll strongly "striped" and anisotropic. Full 
isotropy will be restored only after a second set of simi- 
lar transitions in 0y^ the flrst to floating incommensurate 
order in 0y^ and then a dislocation unbinding transition 
to the liquid phase C. 



As we noted in Section and Fig^ it is possible that 
the low temperature 4x4 commensurate solid can melt 
directly into the high temperature liquid state. In the 
theories T<i> in Eq. (|2.4|l and in Eq. (|3.2|l . such a transi- 
tion appears possible at commensurate densities at which 
T] = T] ~ 0. So this section will neglect the influence of 
77 and perform a full renormalization group analysis of 
the theory J-'g , including both dislocations and the p-fold 
lock-in field h. 

A key property which aids the analysis is a self-duality 
of the action in which the dislocation fugacity, y, and 
the lock-in field h are interchanged. The origin of this 
self-duality is similar to that discussed in Ref. and 
we present a brief derivation in Appendix^ From this 
analysis we obtain the the partition function J-g, includ- 
ing dislocations, obeys 

Zg{Ki,K2, K3,K4, h, y) = Zg{Ki,K2, K3,K4, y, h), 

(5.1) 

where the couplings Ki were defined in Eq. H4.9|l . 

Aided by this duality, we can now immediately deduce 
the full fiow equations for all the elastic constants K^, 
the vortex fugacity y and the field h from the results in 
Section IIV Bl As always, these equations are valid for 
small y and h and are 

di = (2-^^o)/. 
dy P^Tj 

dKi 
~df 



2tt 



p^h^T cos^ </)e^(<^)/^ 
47r2y2 



T 



(Kf+Ki) sin-' (jje 



2 rh^VW/T 



dK2 

di 



d(j> 



p^h^T sin^ 0e^(«/^ 



T 



dK^ 


SiT^y^ 


di 


T 


dKi 


87r2y2 


di 


T 



K2K3 



KIK4 



2-iT 



)cos2 ^e^(^)/^ 



>sin2 ^e^(^)/^ 



(5.2) 



Here Lq and V{(j)) are defined as in Eq. H4.1()|l . but with 
direct couplings: 



Jo 27r 



(5.3) 







10 



We will confine our analysis of Eq. H5.2|l to the phys- 
ically important case of p = 4. For this case, we first 
searched for an intermediate phase with power-law cor- 
relations: such a phase would obtain if there was a set of 
values of Ki^^ for which both y and h flowed to zero. As 
shown in Appendix ^ there is no such regime of param- 
eters. 

However, the flow equations in Eq. H5.2|l do predict 
a direct second-order transition between the 4x4 com- 
mensurate solid A and isotropic liquid C. This transi- 
tion is described by a manifold of fixed points. In the 
6-dimensional space of couplings, there is 2-dimensional 
manifold of fixed points specified by 



y 



K1K2 
T2 



0. 



(5.4) 



In order to describe the flows near this manifold, it is 
convenient to make a change of variables from y, h, Ki, 
and K2 to 



K 



a = y- 
T 



- h 
2 

TT 




(5.5) 



In these variables, the fixed point manifold is described 
by a = if = ~ K4 ~ 0, while the values of A and 
(3 are arbitrary. All physical properties, including the 
exponents at the second-order critical point, will depend 
upon the bare values of A and (3. We expanded Eqs. (|5.2() 
to linear order in a, K, K3, and K^, and after diago- 
nalization of the flow equations, obtained the following 
renormalization group eigenvalues in this 4-dimensional 
subspace: 



I 



1/2N 



1± 1 



V 



7r2/3(\/A + l/\/A) 
47r^/3^ (\/A + 1/\/a)\ 



(5.6) 



where the last eigenvalue is doubly degenerate. It is ev- 
ident that for all A, /3 there are 3 negative eigenvalues 
and 1 positive eigenvalue. This flow therefore describes 
a conventional second-order transition, with the correla- 
tion length exponent v equal to the inverse of the positive 
eigenvalue. 



VI. CONCLUSIONS 



is remains within the realm of possibility that the under- 
lying physics of density wave ordering in all the cuprates 
is the same, and the distinctions between the experiments 
are entirely due to their distinct locations in our phase 
diagram. 

For easy reference, we conclude by listing the various 
routes the 4x4 solid A can melt into the liquid C. The list 
below is not exhaustive, and omits certain possibilities 
involving strong first order transitions between unrelated 
phases. 

1. Phase A undergoes a second order Pokrovsky- 
Talapov^^ (PT) transition to the incommensurate 
striped phase Bg, as described in Section IIII Bl 
This is followed by a first order transition into phase 
B as also described in Section llll Bl and Fig^l Then 
phase B undergoes a dislocation mediated melt- 
ing transition into phase C which is described by 
Eqs. (|4.12() and H4.14I) . This last transition is in a 
universality class which is nearly, but not exactly, 
KT. 

2. PT transition from phase A to phase Bs as in 1. 
Phase B5 then has a KT transition to a striped 
state which has long range order with period 4 
along one direction, and exponentially decaying 
correlations along the other. This striped set melts 
into C after a second set of similar transitions: a PT 
transition into a state with incommensurate quasi- 
long range order in one direction only, and then 
finally a KT transition in C. 

3. First order transition from A to B as described in 
Section IlII Bl and Fig^ Then, a transition from B 
to C as in 2. 

4. At special commensurate densities, there is a di- 
rect, self-dual, second-order transition from A to 
C, described in Section with continuously vary- 
ing exponents. 
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We have shown that thermal fluctuations on 4 x 4 or- 
dered state lead generically to a rather rich phase dia- 
gram as a function of temperature and carrier concentra- 
tion. This phase diagram generically must have regions 
with incommensurate and anisotropic ordering. Thus it 



APPENDIX A: EVALUATION OF U 

This appendix will present the derivation of Lq in 
Eq. (|4.10|) and its application to the flow equations pre- 
sented in Section We begin by writing the integral 
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explicitly as 



dip 
2^ 



(Al) 



K2 cos^ ip + Ki sin^ (p 



K1K2 



dip 





Ki + if 2 



• 2 9 

sm (/3 cos^ p 



(if 2 — ifi) cos2(^ 
Q!+ — cos-^ 2(p a+ — cos^ 



where 



(A2) 



The second term in the last line of Eq. (|AHl is identically 
zero and we are left with 



Lo = {Ki+K2 
= 2 



2ir 



dp 



1 



Q TT a+ — a cos^ 20 

ifl +if2 



a/ Q!+(a+ — a ) 

^i + if2 



iflif2 (ifl + if2)2 - (ifa + i^4)2 



(A3) 



which can be written in terms of our original coupling 
constants Ki as 



4^2 if 1 {Kl -Ki) + K^jKl - Kl) 



(A4) 



This is not to be confused with Lq from Eq. (|5.2() which 
is given by 



ifi + is:2 



(A5) 



Vifiif2 [{Kl + K2Y - {K3 + K^Y] 

For the direct transition from A to C (Section for 
p = 4 we wish to determine if there is any region close 
to the manifold of fixed points where both h and y are 
irrelevant. We do this by expanding Eqs. I|A4|I and HA5|I 
for small K^ and if4 parameterized by 



if3 - hK2 



K4 = 64K1, 



(A6) 



where [Ssl and \Si\ < 1 and substituting into the expres- 
sions for dh/d£ and dy/d£ in Eq. H5.2|l . Define 



it is then straightforward to show that (3{h) ^ if 



if 1 ^ if* 



1 



'■2"3 

7r2if 2 + 4T2 



+ 0(<5f) (A8) 



1.5 



0.5 









— (3(h) = 
- (3(y) = 






P(h) > 

p(y) < 






\ p(h)>6~-, 

\p(y)>o 






P(h)<0 

p(y)>o 









0.5 



le/T 



1.5 



FIG. 5: The signs of P(h) and I3{y) in the ifi — if2 plane for 
if 3 = — 3if2/5 and if 4 = 4ifi/5. The line which indicates 
P{h) = always lies below the line describing I3(y) = and 
consequently h and y are never both irrelevant. 



and f3{y) ^ if 

ifi ^ if* 



7r2 i^l^a + A5iT^ 



TT^Kl 



f2TrK2T{64-S3) 



4T2 
2 



^2i^2 



4T2 



+ 0{Sf) (A9) 



where ifj* = 4T2/(7r2if2). The relative signs of P{h) and 
/3(y) are shown in Fig. Eland from Eqs. I)A8|I and l)A9|l we 
have calculated that the two lines describing the zeros of 
P{h) and (3{y) are always separated by 



32T4if2((54-<53)' 



> 



(AlO) 



(7r2if 2 + 4r2)2 

and thus there are no values of ifi and if2 for small if3 
and if4 where the 4-fold anisotropy and vortices are both 
irrelevant. 



APPENDIX B: DUALITY MAPPING 

This appendix will outline the derivation of Eq. (|5.1|) . 

It is useful to first consider in Eq. (|3.2|l with h = rj = 
0, but to include the effect of vortices. In other words, 
we do wish the impose periodicity under di ^ Oi -\- 27r. 

First, we write the relevant portion of the action in the 
form 

i j d\d,daQa,jbdj9, (Bl) 

where the indices i, j, a, b extend over x, y, and the 
matrix of couplings C can be easily related to the elastic 
constants Ki in Eq. (|3.2I) . Now, we decouple this by a 
set of currents J,v, and write the action as 



(B2) 
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Imposing periodicity in 9i ^ 9i + 27r is now equivalent 
to the requirement that the J^^ are integers. We now 
integrate the 9a out, and solve the resulting constraint 
equations by writing 



2tt 



(B3) 



The constraints that the Jia are integers can now be im- 
posed by demanding that the 9a take values which are 
integer multiples of 2'K/p. As usual, we can soften this 
constraint by introducing a vortex fugacity ?/, and so ob- 



tain the effective action 



ia,jb^3j''^j' 



y^cos{p9a) 



(B4) 

Upon explicitly working out the inverse of the matrix C, 
we find that this action has exactly the same form as the 
original J^g, with the vortex fugacity y playing the role 
of the lock-in field h, and the couplings Ki replaced by 
the couplings Ki in Eq. H4.9|l . 
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